Magnetically driven ferroelectric order in IN13V2O8 
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Abstract 

We show that for N13V2O8 long-range ferroelectric and incommensurate magnetic order appear 
simultaneously in a single phase transition. The temperature and magnetic field dependence of 
the spontaneous polarization show a strong coupling between magnetic and ferroelectric orders. 
We determine the magnetic symmetry of this system by constraining the data to be consistent 
with Landau theory for continuous phase transitions. This phenomenological theory explains our 
observation the spontaneous polarization is restricted to lie along the crystal b axis and predicts 
that the magnitude should be proportional to a magnetic order parameter. 

PACS numbers: 75.80.+q,75.25.+z,75.10.-b 
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The coupling between long-range magnetic and ferroelectric order has been studied since 
the 1960s. 0,0,0] Although a number of systems which are ferroelectric at high temperatures 
become magnetically ordered at a lower temperature, 0, 0]the simultaneous appearance of 
both kinds of order at a single phase transition is much less common. For the most part, 
studies of these multiferroic materials have focussed on commensurate magnets. The magne- 
toelectric properties of these systems are typically discussed in terms of magnetic symmetry 
groups. The onset of magnetic ordering requires breaking time-reversal symmetry, while 
ferroelectric order breaks spatial inversion symmetry. Therefore, only magnetic groups hav- 
ing the proper symmetries allow the possibility of simultaneous magnetic and ferroelectric 
order. 

Recent studies of systems underoing simultaneous ferroelectric order at a magnetic tran- 
sition have identified new multiferroic compounds, including the perovskite Rare Earth 
manganites TbMnC>3and DyMnOq 5], and TblVfr^OsQ. However, the magnetic structures 
of these materials are complex 3, S 2 1 , which makes an investigation of their magnetoelectric 
properties based on symmetry analysis problematic. In contrast, the magnetic structure of 
M3V2O8 (NVO), which we have identified as a multiferroic material, is better determined. 
There are extensive neutron data on this compound |lOj], and a symmetry analysis of the 
data (discussed below) constrains the symmetries which are consistent with a continuous 
transition. Here we take advantage of our analysis of the magnetic symmetry of NVO to de- 
velop a phenomenological Landau theory which provides an explanation of our observations 
of ferroelectric order induced by magnetic ordering. This provides an alternate route to 
understanding the magnetoelectric coupling in multiferroic materials beyond the traditional 
magnetic symmetry group analysis. 

We briefly review earlier results 



10, ystal) of 

(NVO) since its symmetry is crucial to the development of our model. NVO is a magnetic 
insulator consisting of planes of spin-1 Ni 2+ ions arranged in a Kagome staircase lattice. Fig. 
la shows the positions of the two kinds of Ni +2 spins which we call "spine" and "cross tie" 
spins. Competition between several weak magnetic interactions and anisotropics yields the 
complex phase diagram of Fig. 2 and the variety of magnetic structures illustrated in Fig 
lb-d. Cooling at H = 0, one first enters a high temperature incommensurate 3](HTI) phase 
at Tn — 9. IK (Fig. lb), then a low-temperature incommensurate (LTI) phase at = 6.3K 
(Fig. lc). Below 3.9 K the system displays two slightly different canted antiferromagnet 
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(CAF) phases (Fig. Id). The transitions involving the HTI phase are continuous, whereas 
that from the LTI to one of the CAF phases is discontinuous. In the HTI phase the long- 
range order is mostly on the "spine" sites with their spins parallel to a, while in the LTI 
phase the spine and cross tie spins rotate within an a-b plane as shown in Fig. lb-c. 
We used representation analysis [1J] to determine the symmetry of these phases rather than 
attempting to acquire the vast amount of data which would be necessary to experimentally fix 
the several complex order parametersf e. g. see Eq. (2), below] for these magnetic structures. 
Since the irreducible representation associated with magnetic ordering are one dimensional, 
each magnetic order parameter is characterized by a phase factor of unit magnitude which 
specifies how the order parameter transforms under the symmetry operations which leave 
the incommensurate wavevector invariant. Thus the HTI order parameter, o"h is odd under 
a two-fold rotation about the a-axis and is even under the glide plane whose mirror plane 
is an a-b plane. In the LTI phase an additional order parameter, ctl appears which is even 
under both operations. 

Single crystals of NVO were grown from BaO-V 2 5 fluxQ. The single crystal samples 
were oriented using a Laue X-ray diffractometer. To probe the ferroelectric order, we sput- 
tered two gold electrodes on opposite faces of the crystals and measured both the pyroelectric 
current at fixed H, and the magnetoelectric current at fixed T, using a Keithely electrometer. 
We then integrated the current to find the spontaneous polarization P, corresponding to the 
ordered moment arising from the ferroelectric state. We aligned the ferroelectric domains 
by applying a (~ 2 kV/cm) polarizing electric field, E , as the sample was cooled throught 
he transition temperature. After removing the electric field, we measured the polarization 
P a , Pb, and P c , along each of the three crystal axes. To within experimental error, only 
was found to be nonzero. 

Figures 3a and 3b show the electric polarization along the b axis, versus T at zero 
applied voltage for fixed magnetic fields applied along the a and c axes, respectively. Since 
the sign of the polarization in zero electric field changed when the sign of the polarizing 
field-So was changed, we infer that this polarization arises from ferroelectric order in NVO. 
From data for the magnetic field along different crystal axes, it became clear that the region 
in which was nonzero coincided exactly with the region in which the LTI phase existed, 
as shown in Figure 2. Furthermore the hysteresis in ferroelectric order as a function of 
magnetic field or temperature is connected with the fact thatthe LTI-CAF phase transition 
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is discontinuous. 

The isothermal data shown in Fig. 3b and 3d, corroborate the above picture. Fig. 3d 
shows Pf, versus the magnetic field along c, H c , at T = 5K. Atlow H c , Pf, is insensitive to 
the external magnetic field. AsH c is increased, the sample undergoes CAF ordering, which 
completely suppresses the spontaneous polarization. On decreasingif c , returns to the 
initial value. The field-hysteresis is attributed to the first order LTI to CAF transition. Fig. 
3b shows Pb at T =2K, versus H||a. At small H a , there is no ferroelectric order in the CAF 
phase. Increasing H a however produces a spontaneous polarization in the LTI phase which 
is independent of the sign of the magnetic field. The crucial observation is that we can 
magnetically gate NVO to either suppress or promote ferroelectric order, depending only on 
the direction of the applied field and the temperature. 

We now give a phenomenological explanation of our results. In the HTI phase the exis- 
tence of incommensurate magnetic order does not induce a nonzero value of P in the HTI 
phase because the HTI phase is inversion symmetric. This symmetry is not easy to establish 
directly by analyzing the neutron diffraction data. If one assumes the spin amplitudes are 
restricted according to a representation analysis, ^| one still has to determine the complex- 
valued order parameters <JH-i{k) = &H;i{~k)* which characterize the HTI phase, where k 
denotes the wavevector of the incommensurate ordering and i = 1, 6, corresponds to the 
six spin components (the a, b, and c components of the spine and cross tie spins). In the LTI 
phase the additional order parameters (JL;i(k) = crm{—k)* appear. Representation analysis 
is based on the so-called "little group," Gk, which leaves the wavevector invariant. But this 
analysis does not take into account any additional symmetries (such as spatial inversion I) 
which the crystal possesses. We use the inversion symmetry of the crystal lattice to define 
the order parameters so that they obey 

Fcx x , t (k) = ax-M* > (!) 

where X denotes either H or L. The quadratic terms in the Landau free energy which 
describe the appearance of the HTI phase are 

F = Y, F ^H-A-k)a H]j {k) , (2) 

where Fij = F*^ Then Eq. (J!Q) implies that the coefficients F^ are all real valued. Since the 
spin components are those of the critical (i. e.the one whose eigenvalue becomes negative 
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as T is lowered through Th) eigenvectorof the real symmetric matrix F, the amplitudes an,i 
are all real valued, apart from a common overall phase: 



where the q, obtained from the critical eigenvector, are real. Thus the HTI phase is 
described by a single complex valued order parameter e^ H o"H, whose overall phase, 0h, is not 
fixed at this level of analysis. The overall phase (/>h can be eliminated by redefining the origin 
of coordinates. The fact that the resulting Fourier components (TH;i(k) are all real means 
that the magnetic structure is invariant under spatial inversion. Because the relative phases 
of the spin components are difficult to determine from an analysis of the diffraction data, 
this formal argument is essential to establish that the HTI phase is inversion symmetric 
and therefore that the magnetic ordering can not induce a spontaneous electric polarization. 
Since the CAF phaseis also inversion symmetric (Fig. Id), ferroelectricity is not induced 
by magnetic order in that phase either. These results confirm that this model based on our 
extension of representation analysis successfully predicts the experimental observation that 
there is no ferroelectric order in the HTI or CAF phases. We can also understand why when 
a sufficiently large magnetic field is applied so as to enter the CAF phase (Fig. 2b), the 
spontaneous polarization abruptly disappears. 

We now introduce a phenomenological model to explain the symmetry of the magneto- 
electric coupling which induces the spontaneous polarization. In principle, this symmetry 
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might be deduced from a study of magnetic symmetry, jl| but we have found it much simpler 
to directly obtain these results using Landau theory combined with representation analysis. 
While the specific details of this discussion apply to NVO, the general approach is applicable 
to other systems. To describe the incommensurate phases and ferroelectric behavior in a 
single model, we write the Landau free energy as 



where a, and b are constants and xe is the electric susceptibility. (We assume that the 
wavevectors of the two incommensurate phase are locked to be equal.) In a conventional 
ferroelectric P becomes nonzero when \e becomes infinite, often associated with a structural 



<J H -i{k) = ae^an , 



(3) 



F = a(T- T H )4 + b(T - T L )at 




(4) 
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phase transition. Here xe is finite and the appearance of a nonzero P is due to the term 
V coupling magnetic and ferroelectric orders whose form we discuss below. Note that this 
expansion is expressed in the disordered phase, so that all terms in this equation must 
be invariant under the complete set of symmetry operations of the disordered phase. In 
particular, this expression must be invariant under spatial inversion. 

We now discuss the form of the magnetoelectric coupling V. To conserve wavevector 
this coupling must be at least trilinear,|l5( being proportional to one order parameter at 
wavevector k, another at wavevector (—k), and P which is associated with zero wavevector. 
As we have already argued, this coupling is zero with only a single magnetic order parameter 
o"h- A similar argument made for the LTI variables shows that ferroelectricity can not be 
induced if we only invoke the single order parameter ol- Accordingly, and in analogy with 
what has been done for the theory of second harmonic generation, [l6| we posit the following 
interaction which involves two different symmetry order parameters 



V = - ^2{ a iji a H-Ak)a L]j (-k) 



i.r. 



+a* jj cr H ;i(-k)a L; j(k)]P 1 . (5) 

Using I<Jx-i{k) = &x;i{—k) and IP 1 = — P 7 , we see that the inversion invariance of V implies 
that the coefficients must be purely imaginary. Using Eq. (JHJ) and its analog for the L 
order parameters we then have 

Vlti = ^ a -y a HcrL sm((p H - 0l)P 7 , (6) 

7 

where real valued coefficient. At fourth order in the Landau expansion it can be 

shown that tfiu — <pL = tc/2, but to have P 7 7^ it is only essential that 4>h 7^ <Al-(A simple 
way to understand the effect of the fourth order terms is to recall that Wilson in his original 
formulation of the renormalization group, used these terms to mimic the constraint of fixed 
length spins. [12| Here a fixed length of the spins is most nearly achieved by having the two 
incommensurate types of order be out of phase by 7r/2). We now insert Eq.© into Eq. (J3J) 
and minimize with respect to P. We then find a spontaneous polarization given by 

P 7 oc a 7 XeiO"xO"Jf- (7) 

Note that this prediction can be checked by measuring all the quantities which appear in 
Eq. (7) as a function of magnetic field and temperature. This result on the dependence of 



the spontaneous polarization on the magnetic order parameter could not be obtained from 
magnetic symmetry group analysis alone. 

To analyze the consequences of the trilinear coupling in Eq. (JBJ), it is necessary to know 
how the order parameters transform under the symmetry operations of the crystal. As stated 
above, <tl is even and o"h is odd under a two-fold rotation about the a axis, so for V to be 
an invariant, P must also be odd under this operation. (This implies that in Eq. ©, 7 can 
only be b or c.) Furthermore, <tl and <th are both even under the a-b glide plane, which 
restricts 7 in Eq © to be in this plane. Taking account of both symmetries, we see that a 7 
in Eq. (JHJ) can only be nonzero for 7 = 6. This simple argument therefore explains why the 
spontaneous polarization is confined to the b axis. Furthermore, since a) the polarization 
depends on the product of and an and b) an is almost constant in the LTI phase, we 
expect Pb to be approximately proportional to o"l, which is qualitatively what is observed. 
Note also that we do not expect any coupling between the uniform applied magnetic field, 
H a , and the uniform polarization, P 7 , which explains why P 1 was observed to be independent 
of the sign of H a . 

In summary: we have shown that the development of ferroelectric order is coincident 
with an incommensurate magnetic phase in NVO. Since ferroelectricity occurs only in the 
phase for which magnetic ordering breaks inversion symmetry, one can reversibly switch 
the polarization on and off using an external magnetic field. We develop a Landau-like 
model to explain how ferroelectricty is induced by incommensurate magnetic ordering.lt is 
expected that this model will explain the behavior of other recently studied multiferroics with 
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incommensurate or long wavelength structure 18, Ui| and will point the way to microscopic 
theories for these materials. In addition, we suggest that the spontaneous development of 
electric and magnetic long-range order may have been overlooked in simple antiferromagnets 
in which the lattice structure has inversion symmetry but the magnetic moments only have 
inversion symmetry relative to a different origin (Fig le). 
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FIG. 1: Crystal and magnetic structures of NVO. aCrystal structure showing spin- 1 Ni 2 + spine 
sites in red and cross tie sites in blue.b-d Simplified schematic representation of spin arrangement 
in the antiferromagnetic HTI, LTI, and CAF phases. [10] "O" indicates the direction of uniform 
magnetization distributed over spine and cross tie sites in the CAF phase. Only the HTI and CAF 
phases have inversion symmetry relative to the indicated central lattice point. Panel e illustrates 
the center of spatial inversion for the lattice does not coincide with the inversion centre of the 
magnetic structure. 



FIG. 2: Phase diagram of NVO versus Tand H for H||a and H||c in panels a and b respectively. 
The data points indicate anomalies in specific heat (C), magnetization (M), dielectric permittivity 
(e), and electric polarization(P) traces versus H and T. Solid lines are guides to the eye. The 
phases are described in the text and illustrated in Fig. lb-d 



FIG. 3: Promotion and suppression of electric polarization by applying magnetic fields in NVO. 
Temperature and magnetic-field dependence of electric polarization along the b axis for H along 
the a (frames a and b)and c (frames c and d) axes. 
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